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We consider the long-time behaviour of spherically symmetric solutions in the Einstein-Skyrme 
model. Using nonlinear perturbation analysis we obtain the leading order estimation of the tail in 
the topologically trivial sector (B — 0) of the model. We show that solutions starting from small 
compactly supported initial data decay as t~ 4 at future timelike infinity and as u~ 2 at future null 
infinity. We also verified that long-time behaviour for the tail in Einstein-Skyrme model is exactly 
the same as it was obtained for wave maps. 



PACS numbers: 03.50.Kk, 03.65.Pm, ll.10.Lm 



I. INTRODUCTION 

o i 

This paper concerns the late-time asymptotic behaviour of a spherically symmetric self-gravitating Einstein-Skyrme 
(ES) model. It is an extension of the paper Q where we studied quasinormal modes in intermediate asymptotics. It 
is also an extension of work done in Q where the expression for the tail in flat space was obtained. The results of this 
paper are closely connected to the results of paper [5[ where the evolution of wave maps was studied. As we remarked 
in paper , in gravitating Skyrme model the linear perturbation method predicts power-law index 7 = 5 for the tail. 
This estimation is in clear conflict with early numerical results for the tails in ES [l| which were later confirmed by 
the results of paper 2] . To explain this disagreement we have studied the expression for the tail in a gravitating wave 
maps model for details see paper Q where we expected similar long-time asymptotics as for the Skyrme model. In 
the current paper direct calculations in gravitating Skyrme model are performed. 
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In the self-gravitating Skyrme model the most interesting problem is certainly the description of the relaxation 
to the static Skyrmion. Unfortunately, due to the lack of analytic formulae describing static Skyrme soliton, the 
description of this problem is very tedious. To avoid these difficulties we follow in the same way as was done by 
' Bizon et al. [H so that we study the relaxation to the vacuum in the topologically trivial B — sector. To estimate 
the parameters of the tail we apply perturbation techniques elaborated in [Q-Q. Using these techniques we will 
demonstrate that the third-order expression for the tail agrees perfectly with numerical results for small initial data. 
The plan of this paper is as follows. In section II we remind the reader the field equations of the model and shortly 
demonstrate the iterative scheme. We present the difference between Einstein-Skyrme model and wave maps model 
which was analysed in paper pj for t = 1. In the last section we demonstrate the numerical evidence confirming our 
analytical estimations for the tails. 
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II. THEORETICAL BACKGROUND 



■ We consider the Einstein-Skyrme model with dynamics given by the Lagrangian [lfj |: 

L = l^TriVaU^U- 1 ) + J-Tr[(y a U)U-\ (V^tf" 1 ] 2 - — ^fl, (1) 
4 32e^ 167rG 

where V Q is the covariant derivative with respect to the spacetime metric, G - is gravitational constans and R - is a 
scalar of curvature. We assume spherical symmetry and parametrize the metric as follows: 

ds 2 = -e- 2S ^N(r,t)dt 2 + N-^r^dr 2 + rW, (2) 

where dfl 2 is a metric on the unit 2-sphere. Applying the standard hedgehog ansatz U = exp(ia ■ rF(r,t)), where 
a is the vector of Pauli matrices and f - unit radial vector, we obtain the following set of ES equations: 

m = ae- s N 2 PF', (3) 



a ( n ■ 2 t-, sin 4 F ^ T ,P 2 wo. 

- 2 sin 2 F + — + uN{— + F 

2 V r 2 u 2 



(4) 
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(5) 
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P = {e- 5 NuF')' + sm(2F)e- s ( - F' 2 ) - ^4 - 1^ . (6) 

y U r A ) 

Here P and u are auxiliary variables denned as: P = ue s N~ x F and u — r 2 + 2sin 2 i ? , m(t,r) is the mass function 
defined as: m(t, r) = r ( 1 ~ N '> anc [ a = 47rG/ 2 is dimensionless coupling constant. The expression for the tail for a = 
(flat space) was obtained in paper [4|; here we consider gravitating case a > 0. 

To obtain the estimation of the tail we study the evolution of the system described by I©-© starting with small, 
smooth and compactly supported initial data 

F(0,r)=e/(r), F(0, r) = eg(r) . (7) 

Following we postulate perturbation expansion 

m(t,r) = m (t,r) + smi(t,r) + e 2 m 2 (t,r) + . . . , (8) 
S(t,r) = 6 Q (t,r)+e5 1 (t,r)+e 2 S 2 (t,r) + ..., (9) 
F(t,r) = F Q {t,r)+eF 1 (t,r)+e 2 F 2 (t,r)+e 3 F 3 (t,r) + .... (10) 

Collecting the terms with the same power of e we obtain a set of equations for the expansion functions which we solve 
recursively. We are studying the relaxation process to the Minkowski space-time, so too = 8o = Fq = 0. 

In the first order in e the requirement of regularity of the metric function N at the origin and choice of gauge 
8(t, r = 0) = require that mi = 8i = 0. In this perturbation order we obtain free I = 1 radial wave equation for the 
F\ function: 

□f 1 = o, □ = (ii) 

with initial data fi(0, r) = f(r), F\ (0, r) = g(r). The general regular solution of equation (ITT]) has the form 

a '{t-r) + a'{t + r) a(t - r) - a(t + r) 
Fi{t,r) = 1 ^ , (12) 

where the generating function a(r) is determined by initial data. 

In the second perturbation order we obtain the free I = 1 radial wave equation \3F 2 — 0; however, contrary to the 
previous F\(t,r) case, the initial data for F 2 are zero so F 2 has to vanish. In this order of perturbation expansion, 
the metric functions satisfy the following equations 

m> 2 = \r 2 [f 2 + if + ^ , (13) 

m 2 = otr 2 F 1 F[, (14) 
8' 2 = -ar{F 2 +F[ 2 ). (15) 

Finally, in the third order in e we get following equation for F3 

nF Skyrme = ^ J, _ ^ _ ^ _ 2_ ^, p , + + ^2 |^_4 ^ _ 6_ p , _ ^ 



r \ / r \t r 

4 
3r 2 



—r Fi + b , (16) 



where 



b = ^ (*i - 2rFf F{ + r 2 F x (F[ 2 - F 2 )) . (17) 



We may compare this equation for the third-order perturbation F 3 k y rme (t^ r ) w ith corresponding equations obtained 
for wave maps (equation (21) in paper Q) and Skyrmion in flat space (equations (17)-(18) in paper @). To make 
this comparison simpler we rewrite the expression for the wave maps in metric parametrization which we used in this 
paper (see equation ([2])). For 1 = 1 this expression reads: 

nF™ ave ~ map = -25 2 F 1 - <j 2 F x - 5'F[ - - ( m'F! + irnA) + — (4*1 - ~ F l ~ 4i T 

r \ J r \r z r 

+ ^2 F i- (18) 
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We observe that the equation for F^ kyrme (t, r) in gravitating Skyrme model is a generalization of formulae obtained 
for wave maps and flat Skyrme model (see /i-term in equation (18) in paper Both expressions (|16j) and (|18j) are 
linear inhomogeneous wave equations and the only difference is in the form of the source term. Comparing Skyrme 
and wave maps case we observe that the difference appears in the additional inhomogenity in Skyrme case which we 
denoted as b - term. We also immediately see that if we drop b - term in the equation (p~6]> we will get the same 
expression for the tail as it was obtained for the I = 1 wave maps model [Bj]. This expression reads: 



F 3 (t,r) = 



{t 2 - r 2 ) 2 



nf, O [ - 



(19) 



where 



Ci = ^ / (a"( S )) z a(s)ds. 



(20) 



This form of the expression for F^{t,r) 
(i.e. r = const, t — > oo) @: 



leads to the following formula for late-time tail at future timelike infinity 



F 3 (t, 



aC\rt 



(21) 



We are intrested in calculating this correction caused by b - term. As it was remarked by Bizoh et al. , this term is 
of lower order in comparison with other terms which contribute to the formula (|19l) . Therefore we may expect, that 
the contribution from b - term will generate higher order corrections which, in principle, are of the form A/t 1 ', with 
7 > 4. To check this and eventually estimate the power-law index 7 of possible sub-leading contribution we have 
solved numerically the equation = b and we assume the analytical expression for Fi given by this equation (|12l) . 
The results of these calculations are plotted in Fig. 1. 



timelike infinity 
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FIG. 1: The log-log plot of different contributions to the F(t, ro) for fixed ro — 5 as a function of time. We see that F| solution 
decays faster than any power at future timelike infinity; for a reference we have also plotted a power-law with power-law index 
7 = 6. 



In this figure we have plotted two components of full F solution the Fi part resulting from the generating function 



of the form: a(x) = e exp(-x) and the third order correction F3 generated by b - term 



In fact for the F| we have 

plotted two curves corresponding to two different resolutions used in numerical calculations. For both contributions 
we observe a rising part, which depends on the initial data and falling part, which is more universal. As the generating 
function a(x) is effectively the function with compact support, we see that for long times the F\ part of the signal 
vanishes with time ~ exp{—t 2 ). The main contribution to the tail comes from third-order perturbation of F (see 
formula (|16[) ): to make the figure more transparent we do not plot it here. The most interesting curve in this plot 
shows the F| component. It is rising part is also initial data dependent. On its falling part (t > 6) its goes like Fi 
component - decreases faster than any power. 
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For some larger times (t > 7), the shape of the curve F$ seems to be a power-law. However in our opinion this is 
not a real effect, but rather an artifact of our numerical procedure coming from a sort of "ghost potential" (see Q 
for details). There are two reasons which support such hypothesis. First - this part of the F| curve is very steep (for 
comparison see a power-law with the power-law index 7 = 6). It is doubtful if this curve may represent a tail. In 
addition, the shape of the falling part of F% curve depends on the numerical resolution - the better resolution the 
larger part of this curve exhibits the universal behaviour (i.e. decreases faster than any power). This is a typical 
situation in the case of the " ghost potential" . 

Summing this up - the numerical calculations suggest, that the contribution of the b - term falls faster than any 
power - so probably it is equal zero. 

To verify this hypothesis we estimate the correction resulting from the b - term analytically. To do that we 
solve the equation OF^ = b by applying the standard Duhamel formula for solving an inhomogeneous wave equation 
OF = N(t, r) with zero initial data 



t+r-r 



F(t,r) = ^jdT J pP i (fx)N(r,p)dp. 



(22) 



Here Pg(p) are Legendre polynomials, in our model I = 1. Using null coordinates 77 = r — p and £ = t + p and 
denoting the b kernel by K(F) , 



we obtain: 



K(F) 



F J 



2rF 2 F' + r 2 F (F 12 — F 2 ) 



1 

8^ 



t-\-r t — r 

d£ I ^-v)Pi(^K(F 1 (^ V ))dr ) , 



(23) 



(24) 



|t-r| 



where p — (r 2 + (£ — t)(t — r)))/r(£ — rj). We will assume that the initial data Fi(t,r) are compactly supported, i.e. 
they vanish outside a ball of some radius R. Therefore for t > r + R we can drop the advanced part of Fi(t, r). We 
also change the order of integration in (|2"4")l thus we get: 



oc i+r 

F 3 fe (t,r) = i J d v J \t-v)P l { J1 )K(F[ t!t %rj))dZ. 



(25) 



To calculate (|25|) and obtain the estimate of a F| in timelike infinity, we use the following identity (see paper Q): 

t+r 



From (Tl2"|) we have 



(f - r,) n 



(-1) 
(-1) 



,2{n-2f r e+1 (t-r]) n - 



(2£+l)!! [(t - r?) 2 - r 2 ]"- 1 
2(n - 2f r i+1 



F 



l+2-n e+3-n 
2 ' 2 

£ + 3/2 



t-v 



(2£ + l)\\ t l + n 
F 1 {t 1 r) = - 



l + {£ + n)j + 



o'(u) 



F 1 (t,r) = - a"(u 



a(u) 



a'(u) 



F{(t,r) 



a"{u) 



2a' (u) 2a(u) 



r 

Substituting (|2T)) - ([2^)) into (|25p and expanding the function K in inverse powers of p = (£ — T])/2 we get: 

+ 00 t+r 



(26) 

(27) 
(28) 
(29) 



Fl(t,r) = 2 - J d V J 



r/£ 



3a* (r,) 



f5-f (a' 2 ( V )a( V )) 
ar\ 



O 



(30) 
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Performing the inner integral over £ in (|50|) and using the identity (f2l))) we obtain 



-\-oo 



-2'r 



drj 



F(^(»' 3 <'")) + F( 2 |<''»' 3 <"» + f|(»' 2 (")''("» 



(31) 



We have obtained an expression which gives the expansion of the result in inverse powers of t. The factors multiplying 
this inverse powers of t are integrals of total derivatives of expressions which vanish at the integral boundaries, therefore 
they are equal zero. As a result we get the following estimation for F$ at future timelike infinity: 



o 



(32) 



i.e. we have demonstrated that F| does not contains terms of the form A/t 1 with 7 < 7. We would like to stress 
here that in our analysis we considered only two terms to show their cancellation. However, this method may be 
extended to higher orders. To make it possible we have to expand the identity (f2l)]) to higher orders in 1/t. If we do 
that and proceed in the same way we will get the cancelations of higher order coefficients, what is in agreement with 
the hypothesis that F| asymptotically tends to zero. 

In summary - analytical results support the hypothesis stated on the basis of numerical results. Both analytical 
and numerical results are compatible with the fact that the b - term does not contribute to the asymptotic expression 
for the tail. As a results the Skyrme model and £ = 1 wave map model are the examples of models which although 
are different, have the same long time asymptotics. 



III. NUMERICS 



To verify analytical prediction for the tails obtained in the previous section we have performed numerical studies of 
long-time asymptotics in the Einstein-Skyrme model. To do that we have solved numerically the equations ([S])-© with 
initial data described below. For solving evolutional equations we have used method of lines with 5-point, fourth order 
accurate spatial discretization. We have solved the resulting ODE's with fourth order Runge-Kutta method. To solve 
the costraints, i.e. hamiltonian constrain (Q| and slicing condition (|5|) we have also used fourth order Runge-Kutta 
method. Here we need the values of some functions out of the grid - we have obtained them using spline interpolation. 
To ensure regularity at the origin we impose the boundary conditions F(t, r = 0) ~ r and P(t, r = 0) ~ r. To avoid 
the contamination of results by parts of the solution reflected from outer boundary we have used the size of the grid 
big enough, so the solution stops before the reflected signal reaches the observation point. Finally, to suppress the 
accumulation of round-off errors in late times we have used quadrupole precision. In our calculations we have used 
the initial data generated by the function (see (j?l)- (fT2")) ) so we get: a(x) = sexp(~ x 2 ) for different values of e. We 
have started with comparing the behaviour of solutions in the Einstein-Skyrme model and I = 1 wave maps model. 
We have prepared the same initial data and evolved them in both models. The results of these simulations are shown 
in Fig. 2. In the left panel of this figure we plot solutions F(t,ro) for a fixed observation point ro as a function of 
time whereas in the right panel we plot solutions F(u,Vo)/r for fixed value of Do as a function of u. 



timelike infinity null infinity 




FIG. 2: Left panel: The log-log plot of F(t, rn) vs. t for fixed ro = 5. Right panel: The log-log plot of F(u,vo)/r for fixed large 
advanced time vo — t + r — 1000 as the function of retarded time u — t — r. In both model we use a = 0.03 and e = 1.0. 

We observe from above figure differencies at the begining of the evolution and in the intermediate asymptotics but 
for large-times they dissapear. It means that these two models are different but have the same long-time asymptotics. 
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In Fig. 3 we plot F(t, r) in self-gravitating Skyrme model with a — 0.03 for three different values of e. We see that 
on log-log plots the late-time tails are clearly seen as straight lines. 



timelike infinity null infinity 




FIG. 3: Left panel: The log-log plot of F(t,r) for fixed r=5. Right panel: The log-log plot of F(t,r)/r for fixed large advanced 
time v — t + r = 1000 as the function of retarded time u — t — r. In both panels (dotted line) we see that solutions starting 
from small initial data decay as t -4 at future timelike infinity and as u~ 2 at future null infinity 



To obtain the parameters of the tails we should use the following formula: 

F(t,r) = At^ exp(B/t + C/t 2 ) . (33) 

For the comparison with numerical data it is convenient to define the local power index (hearafter LPI) defined as 
follows 

For the assumed form parametrising the tail (|3"3")l we get the following expression for the LPI: 

n(i,r)= 7 +y + (35) 

In Fig. 4 we plot LPI at r — 5 as a function of 1/t. All curves in this figure correspond to small initial data. We 
see that all lines approach the same power-law index 7 = 4 at the future timelike infinity, so numerical data confirm 
analytical prediction for the decay rate see equation (p?Tj) . 
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FIG. 4: The local power index n(t,5) as a function of f/t. 



In Fig. 5 we plot e~ 3 F(t, r) as a function of initial amplitude. According to the analytical prediction the late-time 
behaviour of this quantity does not depend on the magnitude of initial data. We may observe that for not-too-large 
initial data this is really the case. 
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timelike infinity null infinity 




5 25 125 5 25 125 625 

t u 



FIG. 5: Left panel: The log-log plot of e :i F(t,r) vs time for fixed r=5. Right panel: The log-log plot of e z F{t,r)/r for fixed 
large advanced time v = t + r — 1000 as the function of retarded time u — t — r. 

Summary: Using a nonlinear perturbation method, we have demonstrated that leading order formulae for late- 
time tails in ES model are exactly the same as those obtained for wave maps (if we drop this correction b - term). 
In other words, these two models albeit, in principle, different have the same long-time asymptotics and it is seen in 
Fig. 2. We have also checked hypothesis that this correction b - term disapeared faster than any other power so that 
Fg (see equation ([T7])l is probably equal zero (see Fig. 1). We verified numerically the power-law index for the tail 
t~ 4 at future timelike infinity and u~ 2 at future null infinity. We also confirmed that for the Einstein- Skyrme model 
we will get the same analytical results in the leading order for the tails as was obtained by Bizori et al. Q, because 
this b - term does not contribute to the asymptotic expression for the tail. 
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